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ABSTRACT 

The ACME Spectra project provides absolutely calibrated, mostly empirical spectra of exoplanet host stars for use in analysis of the 
stars and their planets. Spectra are obtained from ground-based telescopes and are tied directly to calibrated ground- and space-based 
photometry. The spectra remain only "mostly" empirical because of telluric absorption, but interpolation of stellar models over the 
gaps in wavelength coverage provides continuous stellar spectra. Among other uses, the spectra are suitable for precisely converting 
observed secondary eclipses (occultations) into absolute flux units with minimal recourse to models. In this letter I introduce ACME's 
methods and present a calibrated spectrum of the nearby, super-Earth hosting star 55 Cancri that spans the range from 0.81-5.05 jim. 
This spectrum is well-suited for interpreting near- and thermal-infrared eclipse observations. With this spectrum I show that the 
brightness temperature of the small, low-mass transiting planet 55 Cnc e is 1950^ 2 * K at 4.5 /jm (cooler than previously reported), 
which corresponds to a planetary flux of 0.44^^ mJy. This result suggests the planet has some combination of a nonzero albedo, 
a moderately efficient redistribution of absorbed stellar irradiation, and/or an optically thick atmosphere, but more precise eclipse 
measurements are required to distinguish between these scenarii. 

Key words, infrared: stars — planetary systems — stars: individual (55 Cnc) — stars: fundamental parameters — techniques: 
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1. Introduction 

Measurements of flux densities are the primary diagnostic of 
conditions in stellar and planetary atmospheres, and are used 
to constrain albedos, sizes, temperatures abundances, tempera- 
ture profiles, and even interior conditions of stars and planets. In 
the last decade observations of transiting planets during eclipse 
have begun to reveal the intrinsic emission of these externally ir- 
radiated bodies (e.g., Deming et al., 2005; Charbonneau et al., 
2005). Though initially confined to space-based observatories 
such as Spitzer (and, more recently, Hubble; Swain et al., 2012), 
more favorable targets and improved observing strategies have 
recently enabled a growing number of ground-based measure- 
ments (e.g., Sing & Lopez-Morales, 2009; de Mooij & Snellen, 
2009). 

Because eclipse measurements are inherently relative, pre- 
cise planetary fluxes require stellar fluxes yet more precise. The 
state of the art is now reaching the point at which this condition 
is no longer met; for example, recent measurements of eclipses 
with Spitzer/IRAC have attained relative precisions of 1 - 4% 
(e.g., Agol et al., 2010; Campo et al., 2011), comparable to the 
instrument's 3% absolute calibration accuracy (Spitzer Science 
Center, 2012). 

The goal of the ACME project is to provide a catalog of 
Absolutely Calibrated, Mostly Empirical (ACME) spectra of 
exoplanet host stars with high accuracy and precision. In the 
rapidly approaching era of exoplanet spectroscopy (e.g., Swain 
et al., 2012) it is imperative that such an effort advance beyond 



standard photometric calibrations and focus on calibrated spec- 
troscopy. In addition to the aforementioned conversion of plan- 
etary eclipse depths into absolute flux units, other potential uses 
for such a catalog include the study of the planets' host stars: 
abundance analyses, spectral typing, and refined bolometric lu- 
minosities and effective temperatures are all enabled by such 
data. 

I describe the observation, reduction, and calibration of 
ACME spectra in § 2. As a particular case study I present a cal- 
ibrated spectrum of the nearby star 55 Cancri (55 Cnc) in § 3 1 . 
With this calibrated spectrum I provide an improved measure- 
ment of 55 Cnc e's flux at 4.5 fan in light of recent observations 
(Demory et al., 2012), and I conclude with future prospects. 

2. Observations and Methods 

2.1. Summary of Observations 

Observations thus far have been made solely with SpeX (Rayner 
et al., 2003), a cryogenic, cross-dispersed, near-infrared echelle 
spectrograph at the NASA InfraRed Telescope Facility (IRTF). 
Observations and reduction of SpeX data generally follow the 
steps used in compilation of the IRTF Spectral Library (Cushing 
et al., 2005; Rayner et al., 2009); below I review briefly these 



1 Spectral data are available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http : //cdsweb . 
u-strasbg . fr/cgi-bin/qcat?J/A+A/ 
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steps and highlight the differences between my analysis and that 
of the IRTF Library's. 

All objects are observed using SpeX's short-wavelength 
(SXD) mode. This mode offers continuous wavelength cover- 
age from (except for the gap at 1.81-1. 87 pari) at a 
spectral resolution of 2,000 for the standard 0.3" slit. For bright 
targets such as 55 Cnc (K s = 4), wavelength coverage can be 
extended into the thermal infrared using one of the several long- 
wavelength (LXD) modes, which offer a resolution of 2,500. 

For each target, observations of a nearby AO V star allows 
telluric calibration (Vacca et al., 2003); these stars are typically 
chosen to be located within 15° and 0.1 airmass unit of the pri- 
mary observation. Observations are taken at the parallactic angle 
to avoid spurious spectral slopes induced by differential atmo- 
spheric refraction. I use the standard SpeX calibration macros to 
generate the dark, fiat, and arc frames used to calibrate the data. 
Table 1 shows a log of the observations of 55 Cnc discussed in 
§3. 

2.2. Initial Spectral Reduction 

I use the excellent SpeXTool software package (Cushing et al., 
2004) to reduce, combine, and clean the data, a process that con- 
sists of the following steps 2 . First, optimal extraction converts 
the differences of pairs of nodded frames into uncalibrated spec- 
tra in flux density units. I use a constant background level and 
large apertures (2-3" radii) in the extraction; note that the par- 
ticular choice of reduction parameters is retained in the headers 
of the FITS files generated by SpeXTool. Individual spectra for 
an object are then combined using a robust weighted mean, and 
the mean target spectrum is corrected for telluric absorption us- 
ing the mean AO V spectrum (Vacca et al., 2003). If an object 
was not observed at the parallactic angle, I correct the spectrum 
for losses resulting from chromatic seeing and atmospheric dis- 
persion. Finally, the individual echelle orders are combined by 
applying small (< 2%) offsets to the individual orders. The fi- 
nal result is a spectrum with well-determined uncertainties at all 
wavelengths. In most cases an additional step is the trimming of 
pixels located in regions of extremely low telluric transmission 
or with S/N < 5. 

Note that unlike objects catalogued in the IRTF Spectral 
Library (Cushing et al., 2005; Rayner et al., 2009) the spec- 
tra of ACME targets are not shifted to zero radial velocity; but 
such corrections can easily be made by using the copious body 
of high-precision radial velocity characterization of these sys- 
tems (e.g., Fischer et al., 2008). Future ACME spectra will be 
corrected for interstellar reddening, but owing to its proximity 
55 Cnc requires no such correction (von Braun et al., 2011). 

2.3. Spectral Calibration 

To place an extracted spectrum on an absolute flux scale I fol- 
low Rayner et al. (2009), whose Eq. 1 computes the correction 
factor Cx for each bandpass X by comparing to catalog pho- 
tometry. These factors provide the ratios by which the uncali- 
brated spectrum should be scaled in each bandpass: multiplying 
the spectrum by the weighted average of the Cx calibrates the 
flux scale of the spectrum. Fainter objects observed at < 2.5 fim 
use only 2MASS photometry (Skrutskie et al., 2006) and the 
spectral response curves from Cohen et al. (2003); targets ob- 

2 The SpeXTool procedures used, in order, are: XSPEXT00L, 
XCOMBSPEC, XTELLC0R, XLIGHTLOSS (when appropriate), and 
XMERGEORDERS. 



served at longer wavelengths (such as 55 Cnc) also use WISE 
Wl and/or W2 photometry (Wright et al., 2010) and the spectral 
response curves from Jarrett et al. (201 1). 

Telluric absorption causes several gaps in the observed spec- 
trum: in addition to several smaller regions, the largest gaps 
span 1.81-1.87 fim, and (for longer-wavelength data) from 2.50- 
2.95 fim and 4.18-4.59 fim. These last two gaps intersect the 
WISE Wl and W2 response curves and so must be accounted for 
before flux calibration can take place; this step is also required 
to best constrain bolometric luminosities and for future appli- 
cation to space-based observations not necessarily restricted to 
telluric windows. As described below, to fill these gaps I use 
the BT-Settl library 3 (Allard et al., 2010), which provides high- 
resolution model spectra across a wide range of T e ff, log g and 
[Fe/H]. 

I interpolate the model spectra to obtain a spectrum with 
the stellar parameters appropriate for the target star; for 55 Cnc, 
these are 5196 + 24 K, logg(cgs) = 4.45 ± 0.01, and [Fe/H] = 
0.315 + 0.03 (Valenti & Fischer, 2005; von Braun et al., 2011). 
Convolution of the resultant model with a Gaussian profile 
brings the model to an approximate resolution of 2,000, and 
computation of Cx factors for the model place it on the same flux 
scale as the observed spectrum. At this point the observed and 
model spectra sometimes exhibit mismatched spectral slopes 
at gap edges. To mitigate these discrepancies I measure the 
weighted mean flux density in the 10 nm adjacent to the gap 
boundary for both model and observed spectra; from these mea- 
surements I compute a weighted linear fit that I apply to the 
model spectrum. Thus the details of spectral features in the 
model are retained, while the spectral slope is adjusted to match 
the observations. I re-compute the weighted mean of the Cx fac- 
tors using the continuous, hybrid spectrum; the uncertainty on 
this quantity provides an estimate of the uncertainty in the abso- 
lute flux calibration. Applying this mean factor to the data pro- 
duces an absolutely calibrated, mostly empirical spectrum of the 
target star. 

The situation is slightly more complicated for 55 Cnc, which 
exceeds the nominal WISE saturation limits in the Wl and W2 
bands. For a star of this magnitude, WISE overestimates the flux 
in the Wl and W2 bands by, respectively, 0.1 mag and 0.8 mag 
when compared to fainter stars (Sec. VI.4 of Cutri et al., 2012, 
updated 2011 April 26). I therefore apply these offsets to the 
reported magnitudes before flux-calibrating the spectrum. 

As an additional check, I measure 55 Cnc's flux at 3.6 fim 
and 4.5 fim in archival IRAC photometry (Spitzer Programs 48, 
33970, 40976, 70076, and 80231; P.I.s G. Fazio, K. Luhman, 
M. Marengo, K. Su, and M. Gillon, respectively). I determine 
the stellar flux densities in these bands to be 7.2 + 0.3 Jy and 
4.3 + 0.1 Jy, where the uncertainties represent the spread in val- 
ues across several observing programs. These flux values exceed 
the WISE catalogue fluxes (after applying corrections for the 
slightly different isophotal wavelengths and for partial satura- 
tion, as described above) by 22%. 

Both the IRAC and WISE measurement pairs are internally 
self-consistent, as demonstrated by the fact that both sets yield 

^■(1) equal to 1.05, near the value of unity expected near 
the Rayleigh-Jeans limit. The deciding factor is the comparison 
with the calibrated model spectrum of von Braun et al. (2011) 
described in the following section. Spectral calibration using 
the WISE photometry provides an excellent agreement with von 



3 Available online at http: //phoenix. ens -lyon. fr/ 
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Braun et al.'s model spectrum, but the agreement becomes worse 
by a factor of ~8 when calibrating with the IRAC photometry. 

The 2MASS K s band measurement of 55 Cnc has a qual- 
ity flag of E, indicating a poorly determined flux; therefore I do 
not use this bandpass for flux calibration. In addition, the final 
calibration suggests that the 2MASS J and H measurements un- 
derestimate the stellar flux by roughly 25%. The calibration is 
dominated by the smaller uncertainties of the WISE photometry, 
and I use J, H, Wl, and W2 in the final calibration. 55 Cnc's cal- 
ibrated spectrum is shown in Fig. 1 and is available in electronic 
form. Comparison of this spectrum with the photometrically- 
calibrated Pickles (1998) spectrum of von Braun et al. (2011) 
indicates a good match between the two, as described below. 

3. Results: Stellar and Planetary Fluxes 

3. 1 . Spectrum and Bolometric Flux of 55 Cnc 

The star 55 Cnc hosts five planets (Fischer et al., 2008; 
Dawson & Fabrycky, 2010), including the small transiting planet 
55 Cnc e (Winn et al., 2011; Demory et al., 2011; Gillon et al., 
2012). Because of these planets, and by virtue of its prox- 
imity, the star has been the subject of detailed characteriza- 
tion. Recently von Braun et al. (2011) used new interferomet- 
ric and archival photometric measurements to determine the 
star's bolometric flux (using models from Pickles, 1998) to be 
(1.227 ± 0.0177) x 10-'°Wm" 2 . The calibrated spectrum of 
55 Cnc shown in Fig. 1 emits a flux consistent with that of von 
Braun et al.'s spectrum from 0.8 1-2.49 fxm: the difference be- 
tween the two is only (1.8 + 7.2) x 10" 12 W m -2 , which indepen- 
dently confirms the calibration of both spectra. 

In the IRAC1 and 2 bands the flux densitites of 55 Cnc 
are 6.153 + 0.006 Jy and 3.444 ± 0.006 Jy, respectively, where 
these uncertainties represent only the precision (not the accu- 
racy) of the measurements. The formal uncertainty on these val- 
ues is in fact dominated by the absolute calibration of the ACME 
spectrum: the weighted mean of the Cx factors for 55 Cnc has 
an uncertainty of 12%. Considering this level of accuracy, the 
stellar fluxes are in good agreement with the stellar flux densi- 
ties inferred solely from the calibrated BT-Settl model: 6.03 Jy 
and 3.78 Jy, respectively. This situation will improve for fainter 
stars (which avoid saturation and are measured more precisely 
by 2MASS and WISE), but even here the close correspondence 
(1.4%) between my and von Braun et al.'s in-band fluxes sug- 
gests that the true accuracy of 55 Cnc's spectrum may be better 
than 12%. 

3.2. Brightness Temperature and Flux of 55 Cnc e 

By dint of its high temperature and its host star's apparent bright- 
ness, and despite its small transit depth, 55 Cnc e is the planet 
whose intrinsic emission spectrum is most amenable to charac- 
terization (Winn et al., 201 1; Gillon et al., 2012). At least in the 
near term, this planet will be our sole source of detailed informa- 
tion as to the nature of small, highly irradiated planets and their 
atmospheres. 

The mass, size, and density of 55 Cnc e are all intermedi- 
ate between the Solar System's terrestrial and gas giant planets 
(Dawson & Fabrycky, 2010; Gillon et al., 2012), so its com- 
position cannot be determined solely on the basis of these bulk 
properties (cf. Valencia et al., 2010). However, models of plan- 
etary interiors suggest that 55 Cnc e hosts a substantial envelope 
of volatile materials (Gillon et al., 2012). The planet's mass and 
radius can be fit with an Earthlike core and an envelope of ei- 



ther a H 2 /He (~ 0.03% by mass) or H 2 (~ 20%) (Valencia 
et al., 2010). Models of atmospheric mass loss for CoRoT-7b 
(a smaller, denser, and hotter planet; Leger et al., 2009) reveal 
that it is difficult for such a small, highly irradiated planet to 
retain substantial amounts of H.2/He (Valencia et al., 2010; ?), 
and models of sublimated rock atmospheres suggest that even 
heavier species will not remain aloft beyond the highly irradi- 
ated dayside (Castan & Menou, 2011). Thus it is unlikely that 
55 Cnc e's volatile component is a H.2/He envelope, and it seems 
more likely to be made up of heavier species. 

At present, only observations during transit and/or eclipse 
can constrain the atmospheric properties of such a planet (e.g., 
Miller-Ricci et al., 2009). The first thermal emission from 
55 Cnc e has recently been reported (Demory et al., 2012) via 
detection with Spitzer/IRAC at 4.5 jum of an eclipse whose depth 
is only 131+28 ppm (Demory et al., 2012). This result provides 
a timely opportunity to demonstrate the utility of ACME spec- 
tra: I now demonstrate that the planet's brightness temperature, 
Tb, is cooler than the initially reported value of 2360 ± 300 K 
(Demory et al., 2012). 

The brightness temperature Tb of an eclipsing planet can be 
determined by solving the following relation: 

Fp _ J B A (T B )S x (A)dA /^\ 2 
F* ~ J F*S x (A)dA UJ 

Here F* A is the stellar flux density, Sx(A) is the system response 
in bandpass X, B,\ is the Planck function, and the other sym- 
bols have their usual meanings. In a typical eclipse observation 
the planet/star radius ratio is known to high precision from tran- 
sits and the planet/star flux ratio is measured from eclipses. The 
ACME Spectra project provides F* A , so Tb is the only unknown 
quantity. 

Combining the eclipse depth (Demory et al., 2012), plan- 
etary properties from Gillon et al. (2012), stellar properties 
from Valenti & Fischer (2005) and von Braun et al. (2011), 
the stellar spectrum shown in Fig. 1, and propagating the un- 
certainties in all these quantities using a Monte Carlo approach 
yields T B = 1950+ 2 |$ K (corresponding to 0.44+^ 2 mJy), lower 
than previously reported by 400 K. If 55 Cnc has an exposed, 
rocky component then this new, lower temperature suggests that 
atomic Fe would be more prevalant in the planet's Na-dominated 
atmosphere than SiO, and in general outgassed materials would 
be rather less common than in the hotter case (Castan & Menou, 
2011;?). 

Using a blackbody for F*. instead of the stellar spectrum I 
obtain the published value of T B = 2360^40, a result which em- 
phasizes the point that a star does not emit the same flux den- 
sity as a blackbody at the stellar effective temperature. However, 
note that in converting T B into a flux density I have assumed the 
planet emits like a blackbody in the IRAC2 band; this is proba- 
bly not the case, but it is forgivable considering our current ig- 
norance of the planet's true emission spectrum. 

This new and improved temperature estimate is mildly lower 
(by 1.4cr) than the 2300 K expected in the zero albedo, zero 
redistribution case discussed by Demory et al. (2012). Kepler- 
10b is another small, low-mass, highly irradiated planet whose 
albedo has been inferred from its optical phase curve to be 
0.48 + 0.35 (Batalha et al., 2011; Rouan et al., 2011). It is in- 
triguing to note that such an albedo for 55 Cnc e would make 
the 4.5 ymi brightness temperature quite consistent with the no- 
redistribution case (1930 K). However, this agreement is prob- 
ably coincidental: in contrast to 55 Cnc e, Kepler- 10b contains 
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few or no volatiles (Batalha et al., 201 1) and is presumably of an 
altogether different nature. 

The large uncertainty on the 4.5 prni eclipse depth precludes 
any unambiguous inferences as to the nature of 55 Cnc e's emis- 
sion spectrum. However, if the brightness temperature reported 
here is characteristic of the temperature in the planet's atmo- 
sphere, some combination of the following three scenarios may 
therefore be required. The planet has a nonzero albedo, caused 
either by clouds (Schaefer & Fegley, 2009; Castan & Menou, 
201 1, but see also Leger et al. 201 1) or perhaps by the planet's 
surface; the planet exhibits some redistribution of absorbed stel- 
lar irradiation, suggesting a dense atmosphere (a more tenuous 
atmosphere or a lava ocean cannot transport the requisite heat 
flux; Castan & Menou, 2011; Leger et al., 2011); or the planet 
has an optically thick atmosphere, of which the IRAC2 observa- 
tions probe a particular layer (one not necessarily indicative of 
the planet's surface or equilibrium temperature). Further obser- 
vations at higher precision are of the utmost important if we are 
to distinguish between these scenarii. 

These results demonstrate the utility of calibrated stellar 
spectra for determining the properties of transiting planets and 
their host stars. In particular, the spectrum presented here should 
be of great utility in planning and interpreting future observa- 
tions of 55 Cnc e with Spitzer/IRAC and the Hubble Space 
Telescope. A forthcoming paper will present the first large batch 
of approximately one dozen ACME spectra, several of which 
extend into the thermal infrared. Secondary eclipses have been 
observed for all these targets, and such calibrated spectra should 
provide a useful tools with which to study these systems. 
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Table 1. SpeX Observations of 55 Cnc 



Spectral 


UT Date 


Spectroscopy 


R 


Exp. Time 


AO V 


Airmass 


Seeing 


Type 




Mode 




(sec) 


Standard 


Range 


(arcsec) 


G9IV 


2012 Apr 22 


SXD 


2000 


400 


HD 71906 


1.22-1.25 


1.0 




2012 Apr 22 


LXD2.3 


2500 


66 


HD 71906 


1.14-1.16 


0.75 




Wavelength [fim] 

Fig. 1 Absolutely calibrated, mostly empirical spectrum of 
55 Cnc. The data in black are measured, while those in gray 
are modeled. At bottom are the response curves from 2MASS 
(dotted; Cohen et al., 2003), WISE (dashed; Jarrett et al., 201 1), 
and Spitzer/IRAC (solid; Spitzer Science Center, 2012). 
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